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t^- ■ Abstract 

Vacuum polarization integrals involve the vector spectral functions which can be experimentally 
determined from two sources: (i) e + e~ annihilation cross sections and (££) hadronic r decays. 
Recently results with comparable precision have become available from CMD-2 on one side, and 
ALEPH, CLEO and OPAL on the other. The comparison of the respective spectral functions 
involves a correction from isospin-breaking effects, which is evaluated. After the correction it is 
found that the dominant irir spectral functions do not agree within experimental and theoretical 
uncertainties. Some disagreement is also found for the 47r spectral functions. The consequences 
of these discrepancies for vacuum polarization calculations are presented, with the emphasis 
on the muon anomalous magnetic moment. The work includes a complete re-evaluation of 

* i 

all exclusive cross sections, taking into account the most recent data that became available in 
particular from the Novosibirsk experiments and applying corrections for the missing radiative 
corrections. The values found for the lowest-order hadronic vacuum polarization contributions 
are 

a had,LO = f ( 684 - 7 ± 6 -°ex P ± 3.6 rad ) 1(T 10 [e+e- -based] , 

° M " 1 (709.0 ± 5.1 cxp ± 1.2 rad ± 2.8 SU (2)) 10~ 10 [r-based] , 

where the errors have been separated according to their sources: experimental, missing radiative 
corrections in e + e~ data, and isospin breaking. The Standard Model predictions for the muon 
magnetic anomaly read 

f (11 659 169.3 ± 7.0 had ± 3.5 LBL ± 0.4q ED +ew) 10~ 10 [e+e" -based] , 
° M ~ \ (11 659 193.6 ± 5.9 had ± 3.5 LBL ± 0.4 QED +ew) lO" 10 [r-based] , 

where the errors account for the hadronic, light-by-light scattering and electroweak contributions. 
We observe deviations with the recent BNL measurement at the 3.0 (e + e~) and 0.9 (r) a level, 
when adding experimental and theoretical errors in quadrature. 
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1 Introduction 



Hadronic vacuum polarization in the photon propagator plays an important role in the 
precision tests of the Standard Model. This is the case for the evaluation of the electro- 
magnetic coupling at the Z mass scale, ck(M§), which receives a contribution Aa\ ia d{M^) 
of the order of 2.8 10~ 2 that must be known to an accuracy of better than 1% so that it 
does not limit the accuracy on the indirect determination of the Higgs boson mass from 
the measurement of sin 2 #yK- Another example is provided by the anomalous magnetic 
moment a M = (g^ — 2)/2 of the muon where the hadronic vacuum polarization component 
is the leading contributor to the uncertainty of the theoretical prediction. 

Starting from Refs. [1, 2] there is a long history of calculating the contributions from 
hadronic vacuum polarization in these processes. As they cannot be obtained from first 
principles because of the low energy scale involved, the computation relies on analyticity 
and unitarity so that the relevant integrals can be expressed in terms of an experimentally 
determined spectral function which is proportional to the cross section for e + e _ annihi- 
lation into hadrons. The accuracy of the calculations has therefore followed the progress 
in the quality of the corresponding data [3]. Because the latter was not always suitable, 
it was deemed necessary to resort to other sources of information. One such possibil- 
ity was the use [4] of the vector spectral functions derived from the study of hadronic 
r decays [5] for the energy range less than 1.8 GeV. Another one occurred when it was 
realized in the study of r decays [6] that perturbative QCD could be applied to energy 
scales as low as 1-2 GeV, thus offering a way to replace poor e + e~ data in some energy 
regions by a reliable and precise theoretical prescription [7, 8, 9]. Finally, without any 
further theoretical assumption, it was proposed to use QCD sum rules [10, 11] in order to 
improve the evaluation in energy regions dominated by resonances where one has to rely 
on experimental data. Using these improvements the lowest-order hadronic contribution 
to was found to be [11] 

a had,LO = (692.4 ±6.2) 1(T 10 . (1) 

The complete theoretical prediction includes in addition QED, weak and higher order 
hadronic contributions. 

The anomalous magnetic moment of the muon is experimentally known to very high 
accuracy. Combined with the older less precise results from CERN [12], the measurements 
from the E821 experiment at BNL [13, 14, 15], including the most recent result [16], yield 

a™* = (11 659 203 ±8) 10~ 10 , (2) 

and are aiming at an ultimate precision of 4 10~ 10 in the future. The previous experimen- 
tal result [15] was found to deviate from the theoretical prediction by 2.6 a, but a large 
part of the discrepancy was actually originating from a sign mistake in the calculation of 
the small contribution from the so-called light-by-light (LBL) scattering diagrams [17, 18]. 
The new calculations of the LBL contribution [19, 20, 21] have reduced the discrepancy 
to a nonsignificant 1.6 a level. At any rate it is clear that the presently achieved experi- 
mental accuracy already calls for a more precise evaluation of aJf d ' L °. 
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In this paper we critically review the available experimental input to vacuum polar- 
ization integrals. Such a re-evaluation is necessary because 

• new results have been obtained at Novosibirsk with the CMD-2 detector in the 
region dominated by the p resonance [22] with a much higher precision than before, 
and more accurate R measurements have been performed in Beijing with the BES 
detector in the 2-5 GeV energy range [23]. 

• new preliminary results are available from the final analysis of r decays with ALEPH 
using the full statistics accumulated at LEP1 [24]; also the information from the 
spectral functions measured by CLEO [25, 26] and OPAL [27] was not used previ- 
ously and can be incorporated in the analysis. 

• new results on the evaluation of isospin breaking have been produced [28, 29, 30], 
thus providing a better understanding of this critical area when relating vector r 
and isovector e + e~ spectral functions. 

Since we are mostly dealing with the low energy region, common to both e + e _ and r 
data, and because of the current interest in the muon magnetic moment prompted by the 
new experimental result, the emphasis in this paper is on a^ ad ' LO rather than Aahad(^I)- 
It is true that the presently achieved accuracy on Atthad(^l) is meeting the goals for the 
LEP/SLD/FNAL global electroweak fit. However the situation will change in the long 
run when very precise determinations of sin 2 # w ? as could be available from the beam 
polarization asymmetry at the future Linear Collider, necessitate a significant increase of 
the accuracy on Aa had (M^) [31]. 

Disclaimer: 'theoretical' predictions using vacuum polarization integrals are based 
on experimental data as input. The data incorporated in this analysis are used as quoted 
by their authors. In particular, no attempt has been made to re-evaluate systematic 
uncertainties even if their size was deemed to be questionable in some cases. However, 
whenever significant incompatibilities between experiments occur, we apply an appropri- 
ate rescaling of the combined error. The analysis thus heavily relies on the quality of the 
work performed in the experiments. 



2 Muon Magnetic Anomaly 

It is convenient to separate the Standard Model prediction for the anomalous magnetic 
moment of the muon into its different contributions, 

fl SM = a QED + a had + ^ (3) 

with 

had had,LO _i_ had, HO _i_ had,LBL ( a\ 

where a^ ED = (11 658 470.6 ±0.3) 10~ 10 is the pure electromagnetic contribution (see [32, 
33] and references therein), a^ ad ' LO is the lowest-order contribution from hadronic vacuum 
polarization, a^ ad ' HO = (—10.0 ±0.6) 10~ 10 is the corresponding higher-order part [34, 4], 



2 



and a™ eak = (15.4 ± 0.1 ± 0.2) lO" 10 , where the first error is the hadronic uncertainty 
and the second is due to the Higgs mass range, accounts for corrections due to exchange 
of the weakly interacting bosons up to two loops [35]. For the LBL part we add the 
values for the pion-pole contribution [19, 20, 21] and the other terms [20, 21] to obtain 

a had,LBL = (g Q ± 3 _ 5 ) 1Q -10 

By virtue of the analyticity of the vacuum polarization correlator, the contribution of 
the hadronic vacuum polarization to a M can be calculated via the dispersion integral [36] 

a^o = ^J ds m R[sh (5 ) 

4ml 

where K(s) is the QED kernel [37] , 

K(s) = x 2 + {1 + X)2 { 1 + ^) (Ml+*)-* + y) + jr^y^ lnx ' ^ 

with x = (1 - + and /3 M = (1 - 4mJ/s)V 2 . In Eq. (5), R(s) = R<®(s) 

denotes the ratio of the 'bare' cross section for e + e~ annihilation into hadrons to the 
pointlike muon-pair cross section. The 'bare' cross section is defined as the measured 
cross section, corrected for initial state radiation, electron-vertex loop contributions and 
vacuum polarization effects in the photon propagator (see Section 4 for details). The 
reason for using the 'bare' (i.e. lowest order) cross section is that a full treatment of 
higher orders is anyhow needed at the level of a M , so that the use of 'dressed' cross 
sections would entail the risk of double-counting some of the higher-order contributions. 

The function K(s) decreases monotonically with increasing s. It gives a strong weight 
to the low energy part of the integral (5). About 91% of the total contribution to a^ ad,LO 
is accumulated at center-of-mass energies y/s below 1.8 GeV and 73% of a^ ad,LO is covered 
by the two-pion final state which is dominated by the p(770) resonance. 

3 The Input Data 

3.1 e + e~ Annihilation Data 

The exclusive low energy e + e~ cross sections have been mainly measured by experiments 
running at e + e~ colliders in Novosibirsk and Orsay. Due to the high hadron multiplicity 
at energies above ~ 2.5 GeV, the exclusive measurement of the respective hadronic final 
states is not practicable. Consequently, the experiments at the high energy colliders 
ADONE, SPEAR, DORIS, PETRA, PEP, VEPP-4, CESR and BEPC have measured 
the total inclusive cross section ratio R. 

We give in the following a compilation of the data used in this analysis: 

• The e + e-^ tt+tt- measurements are taken from OLYA [38, 39], TOF [40], CMD [38], 
DM1 [41] and DM2 [42]. 



3 



The most precise data from CMD-2 are now available in their final form [22]. They 
differ from the preliminary ones, released two years ago [43], mostly in the treat- 
ment of radiative corrections. Compared to the preliminary ones, the new results 
are corrected (see Section 4) for leptonic and hadronic vacuum polarization, and 
for photon radiation by the pions (final state radiation - FSR), so that the mea- 
sured final state corresponds to 7r + 7r~ including pion-radiated photons. The various 
changes resulted into a reduction of the cross section by about 1% below the p peak 
and 5% above. The dominant contribution stemmed from vacuum polarization, 
while the (included) FSR correction increased the cross section by about 0.8% in 
the peak region. The overall systematic error of the final data is quoted to be 0.6% 
and is dominated by the uncertainties in the radiative corrections (0.4%). 

We do not use the data from NA7 [44] as they are known to suffer from a systematic 
bias in the energy scale [45]. All the experiments agree with each other within their 
quoted errors, but the high precision claimed by CMD-2 makes this experiment 
unique and consequently not cross-checked by the others at that level. 

The comparison between the cross section results from CMD-2 and from previous 
experiments (corrected for vacuum polarization and FSR, according to the procedure 
discussed in Section 4) is shown in Fig. 1. Note that the errors bars given contain 
both statistical and systematic errors, added in quadrature (this is the case for all 
figures in this paper). The agreement is good within the much larger uncertainties 
(2-10%) quoted by the older experiments. 

• The situation of the data on the oj and resonances has significantly improved 
recently [46, 47, 48]. The numerical procedure for integrating their cross sections is 
described in detail in Section 8.2. 

• The cross sections for e + e - ^ 7r°7 and 777 not originating from the decay of the uj 
and 4> resonances are taken from SND [49] and CMD-2 [50] data in the continuum. 
They include contributions from p — > 7T°7 and p — > 777. 

• The reaction e + e~^ 7r + 7r~7r° is dominated by the u and intermediate resonances 
discussed above. The continuum data are taken from ND [51], DM1 [52], DM2 [53], 
SND [54] and CMD [55]. 

• The e + e~^ 7r+7r-7r 7r° data are available from M3N [56], OLYA [57], ND [51], 
DM2 [58, 59, 60], CMD-2 [61] and SND [62]. It is fair to say that large discrepancies 
are observed between the different results, which are probably related to problems 
in the calculation of the detection efficiency (the cross sections can be seen in Fig. 9 
shown in Section 6.1). The efficiencies are small in general (~ 10 — 30%) and are 
affected by uncertainties in the decay dynamics that is assumed in the Monte Carlo 
simulation. One could expect the more recent experiments (CMD-2 and SND) to 
be more reliable in this context because of specific studies performed in order to 
identify the major decay processes involved. Accordingly we do not include the ND 
data in the analysis. 
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Figure 1: The cross section for e + e~ — > 7r + 7r~ (7) measured by the different experiments. 
The errors bars contain both statistical and systematic errors, added in quadrature. The 
band is the combination of all the measurements used for the numerical integration fol- 
lowing the procedure discussed in Section 7. 



5 



• The reaction e + e~^ loty is mainly reconstructed in the 7r + 7r~7r°7r° final state and is 
thus already accounted for. It was studied by the collaborations ND [51], DM2 [59], 
CMD-2 [61] and SND [62, 63]. We use these cross section measurements to compute 
the contribution corresponding to the u — > 7r°7 decay mode. 

• The e + e~^ 7r + 7r~7r + 7r~ final state was studied by the experiments OLYA [64], 
ND [51], CMD [65], DM1 [66, 67], DM2 [58, 59, 60], CMD-2 [61] and SND [62]. The 
experiments agree reasonably well within their quoted uncertainties (see Fig. 8 in 
Section 6.1). 

• The e + e~^ n + n~ n + n~ n G data are taken from M3N [56] and CMD [65]. It con- 
tains a contribution from the i]7r + 7T~ channel with rj — > 7r + 7r~7r° which has to be 
treated separately because the 77 decay violates isospin. The other five-pion mode 
e + e~^ 7r + 7r~37T° is not measured, but can be accounted for using the isospin rela- 
tion <T 7r + 7r - 37r o = (j 7r + 7r - 7r + 7r - 7r o/2. The relation is used after subtracting the r]n + n^ 
contribution in the 7r + 7r~7r + 7r~7r° rate. Then the i]7r + 7T~ contribution with rj — > 3tt° 
is added to obtain the full 7r + 7r~37r° rate. 

• For the reaction e + e~^ ujty + ty^, measured by the groups DM1 [68], DM2 [53] and 
CMD-2 [69], a contribution is calculated for uo decaying into 7r°7. The dominant 
three-pion decay already appears in the five-pion final state. 

• Similarly, the contribution for e + e~^ u)2n®, with u — > 7r°7, is taken by isospin 
symmetry to be half of e + e~^ ujty + ty^. 

• The process e+e~^ t]tt+tt- was studied by ND [51], DM2 [53] and CMD-2 [69]. We 
subtract from its cross section the contributions which are already counted in the 
7r + 7T~7r + 7r~7r and 7r + 7r~37r° final states. 

• The cross sections of the six-pion final states 37r + 37r~ and 27r + 27r~27r° were measured 
by DM1 [70], CMD [65] and DM2 [71]. For the missing channel vr+TT^Tr one can 
rely on isospin relations in order to estimate its contribution. If only e + e~ data 
are used, the isospin bound [4] is weak, leading to a possibly large contribution 
with an equally large uncertainty. However, some information can be found in the 
isospin-rotated processes 2 t~ — > v T 37r~27r + 7r° and r~ — > v T 27r~7r + 37r°, where the 
hadronic system has been shown [72] to be dominated by u;27r~7r + and cj7r~27r°, once 
the axial-vector r]2n^n + and r]n^27i° contributions [73] are discarded. An isospin 
analysis then reveals the dominance of the up ± 7i T final state. As a consequence 
the 7T + 7r~47r° channel in e + e~ annihilation only receives a very small contribution, 
determined by the 37r + 37r~ cross section. We include a component for u — > 7r°7. 

• The e + e~^ K + K~ and e + e~^ K^R® cross sections above the <f> resonance are taken 
from OLYA [74], DM1 [75], DM2 [76], CMD [77] and CMD-2 [78]. 

• The reactions e+e~^ Kg and e+e~^ K + K~n° were studied by DM1 [79, 80] 
and DM2 [58]. Using isospin symmetry the cross section of the final state K^K^ir 
is obtained from the relation <J K o K o w o = uk+k-w"- 

2 Throughout this paper, charge conjugate states are implied for r decays. 
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• The inclusive reaction e + e — > K$+X was analyzed by DM1 [81]. After subtracting 
from its cross section the separately measured contributions of the final states KgK®, 
Kg K ± tt t and K^K^ , it still includes the modes K$K$Tr + n-, K^Kln+n- and 
K^K ± n T iT . With the assumption that the cross sections for the processes e + e~ — > 
K° K° (tttt) and e + e~ — > K + K~(tt7t)° are equal, one can summarize the total KKmr 
contribution as twice the above corrected K^+X cross section. Implied by the 
assumption made, it is reasonable to quote as the systematic uncertainty one-half of 
the cross section for the channel K + K~it + it~ measured by DM1 [79] and DM2 [82]. 

• Baryon-pair production is included using the cross sections for pp from DM1 [83] 
and DM2 [84], and for nn from FENICE [85]. 

• At energies larger than 2 GeV the total cross section ratio R is measured inclu- 
sively. Data are provided by the experiments 772 [86], MARK I [87], DELCO [88], 
DASP [89], PLUTO [90], LENA [91], Crystal Ball [92, 93], MD-1 [94], CELLO [95], 
JADE [96], MARK- J [97], TASSO [98], CLEO [99], CUSB [100], MAC [101], and 
BES [23]. Due to their weak experimental precision, the data of 772 are not used 
in this analysis. The measurements of the MARK I Collaboration are significantly 
higher than those from more recent and more precise experiments. In addition, the 
QCD prediction of R, which should be reliable in this energy regime, favours lower 
values, in agreement with the other experiments. Consequently the MARK I results 
on R have been discarded. 

Although small, the enhancement of the cross section due to 7 — Z interference is 
corrected for energies above the J/ip mass. We use a factorial ansatz according to 
Ref. [102, 3], yielding a negligible contribution to a^ ad ' L °. 

The R data in the charm region are displayed in Fig. 2. Good agreement is found 
among the experiments. 

• The narrow cc and bb resonances are treated in Section 8.3. 



3.2 Data from Hadronic r Decays 

Data from r decays into two- and four-pion final states t~ — > z/ r 7r~7r°, t~ — > z/ r 7r~37r° 
and T~ -> z/ T 27r-7r + 7T , are available from ALEPH [5], CLEO [25, 26] and OPAL [27]. 
Very recently, preliminary results on the full LEP1 statistics have been presented by 
ALEPH [24]. They agree with the published results, but correspond to a complete re- 
analysis with refined systematic studies allowed by the 2.5 times larger data set. The 
branching fraction B W7T o for the r — > v T tt^tt (7) decay mode is of particular interest since 
it provides the normalization of the corresponding spectral function. The new value [24], 
B nn o = (25.47 ± 0.13) %, turns out to be larger than the previously published one [103] 
based on the 1991-93 LEP1 statistics, (25.30 ± 0.20) %. 

Assuming (for the moment) isospin invariance to hold, the corresponding e + e~ isovec- 
tor cross sections are calculated via the CVC relations 

1 = 1 47rQ;2 /rr\ 

Ve+e-^n+TT- = —Vtt-ttO , (7) 
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Figure 2: R data in the charm region. The band is the combination of all the measurements 
used for the numerical integration following the procedure discussed in Section 7. 
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The r spectral function iy(s) for a given vector hadronic state V is defined by [104] 



mz 



vv[s 



B(r-^u T V-) dN v 

6 \V u d\ 2 Sew B(t~ — > v T e~ v e ) N v ds 



mi 



2s 



mi 



(10) 



where \V u d\ = 0.9748±0.0010 is obtained from averaging 3 the two independent determina- 
tions [105] from nuclear (3 decays and kaon decays (assuming unitarity of the CKM matrix) 
and Sew accounts for electroweak radiative corrections as discussed in Section 5.1. The 
spectral functions are obtained from the corresponding invariant mass distributions, after 
subtracting out the non-r background and the feedthrough from other r decay channels, 
and after a final unfolding from detector effects such as energy and angular resolutions, 
acceptance, calibration and photon identification. 

It is important to note that r decay experiments measure decay rates that include the 
possibility of photon radiation in the decay final state. Depending on the experiment, 
the analysis may (ALEPH) or may not (CLEO) keep events with radiative photons in 

3 Since the two determinations, |T4d|nucicons = 0.9734 ± 0.0008 and |T4d|kaons = 0.9756 ± 0.0006 are 
not consistent, the final error has been enlarged correspondingly. 



8 



the final state, but all experiments rely on the TAUOLA r decay library [106] to com- 
pute their efficiencies. In TAUOLA charged particles are given a probability to produce 
bremsstrahlung using the PHOTOS procedure [107] which is based on the leading loga- 
rithm approximation valid at low photon energy. Thus the measured spectral functions 
correspond to given final states inclusive with respect to radiative photons in the r decay. 

It should be pointed out that the experimental conditions at Z (ALEPH, OPAL) and 
T(4S) (CLEO) energies are very different. On the one hand, at LEP, the r + r~ events 
can be selected with high efficiency (> 90%) and small non-r background (< 1%), thus 
ensuring little bias in the efficiency determination. The situation is not as favorable at 
lower energy: because the dominant hadronic cross section has a smaller particle multi- 
plicity, it is more likely to pollute the r sample and strong cuts must be applied, hence 
resulting in smaller efficiencies. On the other hand, CLEO has an advantage for the 
reconstruction of the decay final state since particles are more separated in space. The 
LEP detectors have to cope with collimated r decay products and the granularity of the 
detectors, particularly the calorimeters, plays a crucial role. One can therefore consider 
ALEPH/OPAL and CLEO data to be approximately uncorrelated as far as experimental 
procedures are concerned. The fact that their respective spectral functions for the 7r~7r° 
and 27T~7r + 7r° modes agree, as demonstrated in Fig. 3 for 7r~7r°, is therefore a valuable 
experimental consistency test. 

4 Radiative Corrections for e + e~ Data 

Radiative corrections applied to the measured e + e~ cross sections are an important step 
in the experimental analyses. They involve the consideration of several physical processes 
and lead to large corrections. We stress again that the evaluation of the integral in 
Eq. (5) requires the use of the 'bare' hadronic cross section, so that the input data must 
be analyzed with care in this respect. 

Several steps are to be considered in the radiative correction procedure: 

• Corrections are applied to the luminosity determination, based on large-angle Bhabha 
scattering and muon-pair production in the low-energy experiments, and small-angle 
Bhabha scattering at high energies. These processes are usually corrected for exter- 
nal radiation, vertex corrections and vacuum polarization from lepton loops. 

• The hadronic cross sections given by the experiments are always corrected for initial 
state radiation and the effect of loops at the electron vertex. 

• The vacuum polarization correction in the photon propagator is a more delicate 
point. The cross sections need to be fully corrected for our use, i.e. 




(11) 
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Figure 3: Comparison between the shapes of the ALEPH, CLEO and OPAL spectral 
functions normalized to the world-average branching ratio B x7T o . The normalization errors, 
correlated between the shapes, are not contained in the error bars. 
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where dressed is the measured cross section already corrected for initial state radia- 
tion, and a(s) is obtained from resummation of the lowest-order evaluation 



= i — x rr- ( 12 ) 

1 - Aa lcp (s) - Aa had (s) 

Whereas Aai ep (s) can be analytically calculated (here given to leading order) 

Aahad(s) is related by analyticity and unitarity to a dispersion integral, akin to (5), 

ie) 



AaUs) = Re Ids' RW . , (14) 

37T 7 s'(s' — S — IE) 

Ami 



which must also be evaluated using input data. Since the hadronic correction in- 
volves the knowledge of R(s) at all energies, including those where the measurements 
are made, the procedure has to be iterative, and requires experimental as well as 
theoretical information over a large energy range. 

This may explain why the vacuum polarization correction is in general not applied 
by the experiments to their published cross sections. Here the main difficulty is even 
to find out whether the correction (and which one? leptonic at least? hadronic?) 
has actually been used, as unfortunately this is almost never clearly stated in the 
publications. The new data from CMD-2 [22] are explicitly corrected for both 
leptonic and hadronic vacuum polarization effects, whereas the preliminary data 
from the same experiment [43] were not. 

In fact, what really matters is the correction to the ratio of the hadronic cross sec- 
tion to the cross section for the process used for the luminosity determination. In 
the simplest case (for example, DM2 for the 7t + tt~ channel) of the normalization 
to the e + e~ — > process, the vacuum polarization effects cancel. However, 

generally the normalization is done with respect to large angle Bhabha scattering 
events or to both Bhabha and In the latter case, Bhabha events dominate 

due to the t-channel contribution. In the 7r + 7r~ mode, all experiments before the 
latest CMD-2 results corrected their measured processes (7r + 7r~, fi + fi~ and e + e~) 
for radiative effects using 0(a 3 ) calculations which took only leptonic vacuum po- 
larization into account [108, 109]. For the other channels, it is harder to find out 
as information about the luminosity determination and the detailed procedure for 
radiative corrections is in general not given in the publications. 

For all e + e~ experimental results, but the newest 7r + 7r~ from CMD-2 and DM2, we 
apply a correction Chvp for the missing hadronic vacuum polarization given by [110] 

_ 1 - 2Aa had (g) 
Chvp " 1 - 2A ahad (t) ' (15) 
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where the correction in the denominator applies to the Bhabha cross section eval- 
uated at a mean value of the squared momentum transfer t, which depends on the 
angular acceptance in each experiment. A 50% uncertainty is assigned to Chvp- 
For the u and <p resonance cross sections, we were informed that the recent CMD-2 
and SND results were not corrected for leptonic vacuum polarization, so in their 
case we applied a full correction taking into account both leptonic and hadronic 
components. 



• In Eq. (5) one must incorporate in R(s) the contributions of all hadronic states 
produced at the energy yt~s. In particular, radiative effects in the hadronic final 
state must be considered, i.e., final states such as V + 7 have to be included. 

Investigating the existing data in this respect is also a difficult task. In the 7r + 7r~ 
data from CMD-2 [22] most additional photons are experimentally rejected to reduce 
backgrounds from other channels and the fraction kept is subtracted using the Monte 
Carlo simulation which includes a model for FSR. Then the full FSR contribution 
is added back as a correction, using an analytical expression computed in scalar 
QED (point-like pions) [111]. As this effect was not included in earlier analyses, we 
applied the same correction to older 7r + 7r~ data. 

In principle one must worry about FSR effects in other channels as well. For the 
inclusive R measurements it is included by definition. When R is evaluated from 
QCD at high energy, the prediction must be corrected for FSR from the quarks, 
but this is a negligible effect for aJf d ' L °. The situation for the exclusive channels 
is less clear because it depends on the experimental cuts and whether or not FSR 
is included in the simulation. Taking as an educated guess the effect in the 7r + 7r~ 
channel, we correspondingly correct the contributions to a^ ad ' LO from all remaining 
exclusive channels by the factor C F sr = (1-004 ± 0.004) nc where n c is the charged 
particle multiplicity in the final state. 



In summary, we correct each e + e~ experimental result, but those from CMD-2 titi, 
by the factor C ra d = ChvpCfsr- As an illustration of the orders of magnitude involved, 
the different corrections in the 7r + 7T - contribution amount to —2.3% for the leptonic 
vacuum polarization, +0.9% for the hadronic vacuum polarization, and +0.9% for the 
FSR correction. The correction to the mr/ee ratio from the missing hadronic vacuum 
polarization is small, typically 0.56%. Both the vacuum polarization and FSR corrections 
apply only to experiments other than CMD-2, therefore the overall correction to the 7T7T 
channel is considerably reduced. 

The uncertainties on the missing vacuum polarization (50%) and the FSR corrections 
(100%) are conservatively considered to be fully correlated between all channels to which 
the correction applies. The total error from these missing radiative corrections, taken as 
the quadratic sum of the two contributions, is given separately for the final results. 
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5 Isospin Breaking in e + e and r Spectral Functions 

5.1 Sources of Isospin Symmetry Breaking 

The relationships (7), (8) and (9) between e + e~ and r spectral functions only hold in the 
limit of exact isospin invariance. This is the Conserved Vector Current (CVC) property 
of weak decays. It follows from the factorization of strong interaction physics as produced 
through the 7 and W propagators out of the QCD vacuum. However, we know that we 
must expect symmetry breaking at some level from electromagnetic effects and even in 
QCD because of the up and down quark mass splitting. Since the normalization of the 
r spectral functions is experimentally known at the 0.5% level, it is clear that isospin- 
breaking effects must be carefully examined if one wants this precision to be maintained 
in the vacuum polarization integrals. Various identified sources of isospin breaking are 
considered in this section and discussed in turn. 

Because of the dominance of the 7nr contribution in the energy range of interest for r 
data, we discuss mainly this channel, following our earlier analysis [4]. The corrections on 
a had,LO f rom isospin breaking are given in Table 1. A more complete evaluation is given 
in the next section. Finally, the 4-pion modes will be briefly discussed. 

• Electroweak radiative corrections must be taken into account. Their dominant con- 
tribution comes from the short distance correction to the effective four-fermion cou- 
pling t~ — > v T [du)~ enhancing the r amplitude by the factor (1 + 3a(m T )/47r)(l + 
2Q) In (M z /m T ), where Q is the average charge of the final state partons [112]. While 
this correction vanishes for leptonic decays, it contributes for quarks. All higher- 
order logarithms can be resummed using the renormalization group [112, 113], and 
the short distance correction can be absorbed into an overall multiplicative elec- 
troweak correction S^, 

, M _ ( V"» («(M W ) f X ( a(M z) f 17 
bEV >-{ a (m T )) {a(m b )j \a(M w )) ' (lb> 

which is equal to 1.0194 when using the current fermion and boson masses and for 
consistency [114] the quark-level MS expressions for a(s) as given in Ref. [115]. The 
difference between the resummed value and the lowest-order estimate (1.0188) can 
be taken as a conservative estimate of the uncertainty. QCD corrections to 
have been calculated [112, 116] and found to be small, reducing its value to 1.0189. 

Subleading non-logarithmic short distance corrections have been calculated to order 
0(a) at the quark level [117], S s E u w Md = 1 + a(m T )(85/12 - n 2 )/(2n) ~ 0.9967, 
and for the leptonic width [112], S s E n w Acp = 1 + a(m T )(25/4 - 7r 2 )/(27r) ~ 0.9957. 
Summing up all the short distance corrections, one obtains the value for Sew that 
must be used for the inclusive hadronic width 

Chad r<sub,had 
^mdusive = EW^EW = imQg ± Q QQQg _ 

Sew 

Other uncertainties on the b quark mass, the running of a(s), and QCD corrections 
are at the 10~ 4 level. 
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Long distance corrections are expected to be final-state dependent in general. They 
have been computed for the r~ — > u t tt~ decay leading to a total radiative correction 
of 2.03% [118], which is dominated by the leading logarithm from the short distance 
contribution. Although very encouraging, this result may not apply to all hadronic 
r decays, in particular for the important z/ r 7r~7r° mode. Therefore an uncertainty 
of 0.0040 was previously assigned to Sew ( see the following Section 5.2 for more) to 
cover the final-state dependence of the correction with respect to the calculation at 
the quark level. 

A contribution [28, 4] for isospin breaking occurs because of the mass difference 
between charged and neutral pions, which is essentially of electromagnetic origin. 
The spectral function has a kinematic factor f3 3 which is different in e + e~ (it + it~) 
and r decay (7r~7r°). We write 

Ms) = ^\F°(s)\\ (18) 

v-(s) = ^l^)' 2 ' ( 19 ) 

with obvious notations, F®~(s) being the electromagnetic and weak pion form fac- 
tors, respectively, and f3o- defined by 

A) _ = P(s, m w - , 771*0,- ) , (20) 

where 



/3(s,m 1 ,m 2 ) 



(mi + m 2 ) 2 \ ( (mi — m 2 ^ 2 



1-^- ^ 1- 



1/2 



(21) 



Hence, a correction equal to (3q(s) / /3^_(s) is applied to the r spectral function. 

Other corrections occur in the form factor itself. It turns out that it is affected by 
the pion mass difference because the same ft 3 factor enters in the p — > ixix width. 
This effect partially compensates the (3 3 corrections (18), (19) of the cross section, 
as seen in Table 1. 

Similarly a possible mass difference between the charged and neutral p meson affects 
the value of the corresponding width and shifts the resonance lineshape. Theoret- 
ical estimates [119] and experimental determinations [5, 120] show that the mass 
difference is compatible with zero within about 1 MeV. 

p — uo interference occurs in the 7i + n~ mode and thus represents an obvious source 
of isospin symmetry breaking. Its contribution can be readily introduced into the 
t spectral function using the parameters determined in the CMD-2 fit [22]. The 
integral over the interference almost vanishes by itself since it changes sign at the u 
mass, however the s-dependent integration kernel produces a net effect (Table 1). 

Electromagnetic p decays explicitly break SU(2) symmetry. This is the case for 
the decays p — > mv°j through an uj intermediate state because of identical 7r°'s, 
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Sources of Isospin Afl had,Lo ( 10 -io) 



Symmetry Breaking it + it (I) it + it (II) it + it 2n° 2ir + 2ir 

Short distance rad. corr. _ 1Q 3 ± ^ -12.1 ± 0.3 _ Q ^ ± Q Q? _ Q lg ± ^ 

Long distance rad. corr. —1.0 

m?r _ ^ 77V, (/? in cross section) -7.0 -7.0 +0.6 ± 0.6 -0.4 ± 0.4 

m w - ^ r?v> (P in p width) +4.2 +4.2 

m p - ^ m p0 + 0.2 ± 2.0 

p — iv interference +3.5 ± 0.6 +3.5 ± 0.6 

Electromagnetic decay modes —1.4 + 1.2 —1.4 + 1.2 

Sum -11.0 + 2.5 -13.8 + 2.4 +0.2 + 0.6 -0.6 + 0.4 



Table 1: Expected sources of isospin symmetry breaking between e + e and r spectral 
functions in the 2ir and An channels, and the corresponding corrections to 

a had,LO as 

obtained from r data. The corrections (I) follow essentially the procedure used in Refs. [4, 
7, 11], while in (II) the more complete approach of Ref. [30] is chosen. The values given 
for (II) differ slightly from those quoted in Ref. [30], because of the model used in the latter 
to parametrize the pion form factor, in addition to the re- evaluation of the short distance 
electroweak correction. The errors given are theoretical only. Uncertainties introduced by 
the experimental error on the r spectral function itself are not accounted for here. 

p — > 7T7, p° — > 777 and p° — > . The decay p — > 7r7r7 deserves particular attention: 
calculations have been done with an effective model [121] for both charged and 
neutral p's. The different contributions are listed in Table 1. 

• A breakdown of CVC is due to quark mass effects: m u different from rrid gen- 
erates <9 M ~ (m u — md) for a charge- changing hadronic current between u 
and d quarks. Expected deviations from CVC due to so-called second class cur- 
rents such as, e.g., the decay t~ — > z/ T 7r~r/ where the corresponding e + e~ final 
state n°rj (C=+l) is forbidden, lead to an estimated branching fraction of the order 
of (m u — md) 2 jrr? T — 10 -5 [122], while the experimental upper limit amounts to 
B(t -> v T n-ri) < 1.4 10~ 4 [105]. 



5.2 A More Elaborate Treatment of Isospin Breaking in the 2tt Channel 

The above analysis of isospin breaking leaves out the possibility of sizeable contributions 
from virtual loops. This problem was studied recently [29] within a model based on 
Chiral Perturbation Theory. In this way the correct low-energy hadronic structure is 
implemented and a consistent framework can be set up to calculate electroweak and 
strong processes, such as the radiative corrections in the r — > z/ r 7r~7r° decay. One might 
worry that the p mass is too large for such a low-energy approach. However a reasonable 
matching with the resonance region [123] and even beyond is claimed to be achieved, 
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providing a very useful tool to study radiative decays. 

A new analysis has been issued [30] which is more suited to our purpose, in the sense 
that it applies to the inclusive radiative rate, r — > i/ T 7r - 7r° (7), as measured by the 
experiments. A consistent calculation of radiative corrections is presented including real 
photon emission and the effect of virtual loops. All the contributions listed in the previous 
section are included and the isospin-breaking contributions in the pion form factor are 
now more complete. Following Ref. [30], the relation between the Born level e + e~ spectral 
function and the r spectral function (19) reads 



1 ] G EM (s) (3 3 - 



2 



u w - ff o (7 )(s) , (22) 



where (jem(s) is the long-distance radiative correction involving both real photon emission 
and virtual loops (the infrared divergence cancels in the sum). Note that the short- 
distance Sew correction, discussed above, is already applied in the definition of v-(s) (cf. 
Eq. (10)), but its value differs from Eq. (17) because subleading quark-level and hadron- 
level contributions should not be added, as double counting would occur. The correct 
expression for the nn mode therefore reads 

Sew(s) = ™ ™p ' = (1-0233 ± 0.0006) ■ Gem(s) , (23) 

the subleading hadronic corrections being now incorporated in the mass-dependent GW(s) 
factor. The form factor correction is dominated by the effect of the pion mass difference 
in the p width, but it also includes a small contribution at the 10 -3 level from the 'chiral' 
form used for the p lineshape. In practice, however, the correction is independent of the 
chosen parametrization of the form factor. The different contributions to the isospin- 
breaking corrections are shown in the second column of Table 1. The values slightly differ 
from those given in Ref. [30] because the authors use a model for the pion form factor 
rather than integrating experimental data. The largest difference however stems from our 
re-evaluation of the short- distance electroweak correction, .Sew, including the subleading 
leptonic contribution. The sum amounts to (—13.8 ± 2.4) 10~ 10 to be compared with 
(-12.0 ± 2.6) 10~ 10 given in Ref. [30]. 

The dominant uncertainty in this method stems from the p ± -p° mass difference. In- 
deed, in the chiral model used in Ref. [30] the only parameter entering the pion form 
factor is the p mass, since the width is given by V p (s) = m p s f3 3 (s)/(9Girfj;). In the 
method previously used and recalled in Section 5.1, the width at the pole was taken as an 
independent parameter with T p (s) = T p ^/s (3 3 (s) /(rn p (3 3 (m p v )), so that the effect of the p 
mass difference approximately cancels after integration. This explains the large difference 
in the uncertainties quoted for the two evaluations in Table 1. 

Since the integral (5) requires as input the e + e~ spectral function including FSR pho- 
ton emission, a final correction is necessary. It is identical to that applied in the CMD-2 
analysis [22, 111] (cf. Section 4). All the corrections are drawn versus s in Fig. 4. The 
overall correction reduces the r rate below the p peak, but, somewhat unexpectedly, has 
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Figure 4: Mass-squared-dependent corrections applied to the it it spectral function from 
r data, following the analysis of Ref. [30]. 



the opposite effect above. This behavior is driven by the long-distance radiative correc- 
tions contained in Gem(s). 

The total correction to the r result in this method, not including the FSR contri- 
bution, amounts to Aa]} ad = (-13.8 ± 2.4) 1(T 10 , where the main contribution to the 
error is due to the experimental limits on the p mass difference. After including the FSR 
contribution, it becomes (—9.3 ± 2.4) 10 -10 , a value consistent with the result in the first 
column of Table 1 which does not include the virtual corrections and uses a less sophis- 
ticated treatment of radiative decays. In the following we apply the correction functions 
from the more complete analysis (method (II) in Table 1) and keep the corresponding 
uncertainty separate from the purely experimental errors. 



5.3 Isospin Breaking in 4tt Channels 

There exists no comparable study of isospin breaking in the 4tt channels. Only kinematic 
corrections resulting from the pion mass difference have been considered so far [28] , which 
we have applied in this analysis. It creates shifts of —0.7 10~ 10 (—3.8%) and +0.1 10~ 10 
(+1.1%) for 27r + 27r~ and 7i + 7i^2tt°, respectively. However, since the four-pion contribution 
to 

a had,LO j g re l a tively less important than the two-pion part (by a little more than an order 
of magnitude in the integration range up to 1.8 GeV) and the experimental uncertainties 
are much larger, we feel this is a justified procedure at the present level of accuracy of 
the data. Moreover, the entire correction has been attributed as systematic error which 
is kept separate from the experimental errors on a^ ad ' LO from these channels. 
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It should also be pointed out that the systematic uncertainties from isospin breaking 
are essentially uncorrelated between the 2ir and 4ir modes: as Table 1 shows, the dominant 
sources of uncertainties are the p ± -p° mass difference for 2n and the threshold factors in 
47r where large errors have been given to cover uncertainties in the decay dynamics and 
the missing pieces. 

6 Comparison of e + e~ and r Spectral Functions 

The e + e~ and the isospin-breaking corrected r spectral functions can be directly compared 
for the dominant 2-k and 4n final states. For the 2n channel, the p-dominated form factor 
falls off very rapidly at high energy so that the comparison can be performed in practice 
over the full energy range of interest. The situation is different for the 4-7T channels where 
the r decay kinematics limits the exercise to energies less than ~ 1.6 GeV, with only 
limited statistics beyond. 

6.1 Direct Comparison 

Fig. 5 shows the comparison for the 2-7T spectral functions. Visually, the agreement seems 
satisfactory, however the large dynamical range involved does not permit an accurate test. 
To do so, the e + e~ data are plotted as a point-by-point ratio to the r spectral function in 
Fig. 6, and enlarged in Fig. 7, to better emphasize the region of the p peak 4 . The e + e _ 
data are significantly lower by 2-3% below the peak, the discrepancy increasing to about 
10% in the 0.9-1.0 GeV region. 

The comparison for the 4tt cross sections is given in Fig. 8 for the 27r + 27r~ channel and 
in Fig. 9 for 7t + 7t~2tt . As noted before, the latter suffers from large differences between 
the results from the different e + e~ experiments. The r data, combining two measured 
spectral functions according to Eq. (9) and corrected for isospin breaking as discussed in 
Section 5, lie somewhat in between with large uncertainties above 1.4 GeV because of the 
lack of statistics and a large feedthrough background in the r — > v T n~3ir° mode. In spite 
of these difficulties the 7r~37r° spectral function is in agreement with e + e~ data as can be 
seen in Fig. 8. It is clear that intrinsic discrepancies exist among the e + e~ experiments 
and that a quantitative test of CVC in the 7r + 7r~27r° channel is premature. 

6.2 Branching Ratios in r Decays and CVC 

A convenient way to assess the compatibility between e + e~ and r spectral functions pro- 
ceeds with the evaluation of r decay fractions using the relevant e + e~ spectral functions 
as input. All the isospin-breaking corrections detailed in Section 5.2 are included. The 
advantage of this procedure is to allow a quantitative comparison using a single number. 
The weighting of the spectral function is however different from the vacuum polarization 
kernels. Using the branching fraction B(t~ — > v T e~ 9 e ) = (17.810 ± 0.039)%, obtained 

4 Thc central bands in Figs. 6 and 7 give the quadratic sum of the statistical and systematic errors of the 
combined r spectral functions. Local bumps in these bands stem from increased errors when combining 
different experiments having local inconsistencies. We use the procedure described in Section 7.1. 
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Figure 5: Comparison of the 7r + 7r~ spectral functions from e + e~ and isospin-breaking 
corrected r data, expressed as e + e~ cross sections. The band indicates the combined e + e~ 
and t result within la errors. It is given for illustration purpose only. 
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Figure 6: Relative comparison of the 7t + tt~ spectral functions from e + e~ and isospin- 
breaking corrected r data, expressed as a ratio to the r spectral function. The band shows 
the uncertainty on the latter. 




Figure 7: Relative comparison in the p region of the 7r + 7r~ spectral functions from e + e~ 
and is o spin-breaking corrected r data, expressed as a ratio to the r spectral function. The 
band shows the uncertainty on the latter. 
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Mode 


r data 


Branching 

e+e' 


; fractions (in %) 
- via CVC 


A(r - e+e-) 


T 
1 


± 1J TT 7T^ 

' iff i\ J I 






+ n ii ,4-01 2ott^ 


ii 48 + n 32 


T~ 


— > ^ T 7r~37r° 


1.01 ±0.08 


1.09±0.06 eX p 


0.30 

± 0.02 rad ± 0.05 SU (2) 

v ' 


-0.08 ±0.11 




— > ^ T 27r _ 7T + 7T 


4.54 ±0.13 


3.63±0.19 cxp 


0.08 

± 0.04 rad ± 0.09 su(2) 
, ' 


±0.91 ±0.25 










0.21 





Table 2: Branching fractions of r vector decays into 2 and 4 pions in the final state. 
Second column: world average. Third column: inferred from e+e" spectral functions 
using the isospin relations (7-9) and correcting for isospin breaking. The experimental 
error of the 7r + 7r~ CVC value contains an absolute procedural integration error of 0.08%. 
Experimental errors, including uncertainties on the integration procedure, and theoretical 
( missing radiative corrections for e + e~ , and isospin- breaking corrections and V u d for r ) are 
shown separately. Right column: differences between the direct measurements in r decays 
and the CVC evaluations, where the separate errors have been added in quadrature. 

assuming leptonic universality in the charged weak current [24], the results for the main 
channels are given in Table 2. The errors quoted for the CVC values are split into un- 
certainties from (i) the experimental input (the e + e~ annihilation cross sections) and the 
numerical integration procedure, (ii) the missing radiative corrections applied to the rele- 
vant e + e~ data, and (Hi) the isospin-breaking corrections when relating r and e + e~ spec- 
tral functions. The values for the r branching ratios involve measurements [24, 124, 125] 
given without charged hadron identification, i.e. for the /i7r°z/ T , /i37t°z/ t and 3/i7r°z/ T final 
states. The corresponding channels with charged kaons have been measured [126, 127] 
and their contributions can be subtracted out in order to obtain the pure pionic modes. 

As expected from the preceding discussion, a large discrepancy is observed for the r — > 
u T 7r _ 7r° branching ratio, with a difference of (— 1.48±0.12 T ±0.25 ee ±0.11 ra d±0.12gu(2))%, 
where the uncertainties are from the r branching ratio, e + e~ cross sections, e + e~ missing 
radiative corrections and isospin-breaking corrections (including the uncertainty on V U( i), 
respectively. Adding all errors in quadrature, the effect represents a 4.6 a discrepancy. 
Since the disagreement between e + e~ and r spectral functions is more pronounced at 
energies above 750 MeV, we expect a smaller discrepancy in the calculation of a^ ad > LO be- 
cause of the steeply falling kernel K(s) in this case. More information on the comparison 
is displayed in Fig. 10 where it is clear that ALEPH, CLEO and OPAL all separately, but 
with different significance, disagree with the e + e~-based CVC result. 

The situation in the 4it channels is different. Agreement is observed for the 7r~37r° 
mode within an accuracy of 11%, however the comparison is not satisfactory for the 
27T~7r + 7r° mode. In the latter case, the relative difference is very large, (22 ± 6)%, com- 



22 





i 






1 i i i i | i , 

CLEO 






1 


-•— 


1 25.42 ±0.42 










OPAL 






I 


-•- 


— 1 25.44 ±0.34 








(-•- 


A 1 PP — 1 n ro 1 i m i n a r\i 

1 25.47 ±0.1 3 






cvc 




Average 




l • 

, i , 


23.98 ±0.30 

, , 1 




25.46 ±0.1 2 

, I I I I I I I 


23 


24 


25 




26 27 



B(x"->v T JtV) (in%) 

Figure 10: JTie measured branching ratios for r — > v t 7t~7t° compared to the prediction from 
the e + e~ — > 7r + 7r~ spectral function applying the is o spin-breaking correction factors dis- 
cussed in Section 5.2. The measured branching ratios are from ALEPH [24], CLEO [124] 
and OPAL [125]. The OPAL result was obtained from their hn branching ratio, reduced 
by the small Kir contribution measured by ALEPH [126] and CLEO [127]. 



pared to any reasonable level of isospin symmetry breaking. As such, it rather points to 
experimental problems that have to be investigated. 



7 The Integration Procedure 

The information used for the evaluation of the integral (5) comes mainly from direct mea- 
surements of the cross sections in e + e~ annihilation and via CVC from r spectral func- 
tions. In general, the integrals themselves are evaluated using the trapezoidal rule, i.e., 
combining adjacent measurement points by linear interpolation. Even if this method is 
straightforward and free from theoretical assumptions (other than CVC in the r case) , its 
numerical calculation requires special care. The finite and variable distance between adja- 
cent measurements creates systematic uncertainties that have to be estimated. The com- 
bination of measurements from different experiments taking into account correlations — 
both within each data set and between different experiments — is the subject of additional 
discussions presented in the following. 
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7.1 Averaging Data from Different Experiments 



To exploit the maximum information from the available data, we combine weighted mea- 
surements of different experiments at a given energy instead of calculating separately the 
integrals for every experiment and finally averaging them. 

The solution of the averaging problem is found by minimizing 

x 2 = E E « - mr 1 K - kj) , (24) 

n=l i,j=l 

where xf is the ith cross section measurement of the nth experiment in a given final state, 
Cfj is the covariance between the ith and the jth measurement and fcj is the unknown 
distribution to be determined. The covariance matrix C n is given by 

rm _ J (A^stat) 2 + (^sys) 2 for 2 = j „• i AT fo^ 

° y " I A^ ys -A- sys for i / , ' V ' • ( - ,] 

where A" stat (A" sys ) denotes the statistical (systematic) error of xf. The systematic errors 
of the e + e~ annihilation measurements are essentially due to luminosity and efficiency 
uncertainties. It is conservative to take them as common errors of all data points of a 
given experiment. The minimum condition dy 2 jdki = 0, Vi leads to the system of linear 
equations 

EE(^-y^)" 1 = °> < = l,...,iVn. (26) 

n=l j=l 

The inverse covariance C^ 1 between the solutions fcj, kj is the sum of the inverse co vari- 
ances of each experiment 

Nexp 

^ = E (CS-) -1 • (27) 

n=l 

If different measurements at a given energy show inconsistencies, i.e., their x 2 per number 
of degrees of freedom (DF) is larger than one, we rescale the error of their weighted average 

by yJ X 2 /T>F. 



7.2 Correlations between Experiments 

Eq. (27) provides the covariance matrix needed for the error propagation when calculating 
the integrals over the solutions fc, from Eq. (26). Up to this point, CV,- only contains corre- 
lations between the systematic uncertainties within the same experiment. However, due 
to commonly used simulation techniques for the acceptance and luminosity determina- 
tion as well as state-of-the-art calculations of radiative corrections, systematic correlations 
from one experiment to another occur. It is obviously a difficult task to reasonably esti- 
mate the amount of such correlations as they depend on the reconstruction capabilities 
of the experiments and the theoretical understanding of the underlying decay dynamics. 
In general, one can state that in older experiments, where only parts of the total solid 
angle were covered by the detector acceptance, individual experimental limitations should 
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dominate the systematic uncertainties. Potentially common systematics, such as radia- 
tive corrections or efficiency, acceptance and luminosity calculations based on the Monte 
Carlo simulation, play only minor roles. The correlations between systematic errors below 
2 GeV energy are therefore estimated to be between 10% and 30%, with the exception of 
the 7r + 7i~ final state, where we impose a 40% correlation due to the simpler experimental 
situation and the better knowledge of the dynamics which leads to non-negligible system- 
atic contributions from the uncertainties of the radiative corrections. At energies above 

2 GeV the experiments measured the total inclusive cross section ratio R. Between 2 and 

3 GeV, individual technical problems dominate the systematic uncertainties. At higher 
energies, new experiments provide nearly full geometrical acceptance which decreases the 
uncertainty of efficiency estimations. Radiative corrections as well as theoretical errors of 
the luminosity determination give important contributions to the final systematic errors 
quoted by the experiments. We therefore estimate the correlations between the system- 
atic errors of the experiments to be negligible between 2 GeV and 3 GeV, 20% between 
3 GeV and 10 GeV. These correlation coefficients are added to all those entries of CV,- 
from Eq. (27) which involve two different experiments. 



7.3 Evaluation of the Integral 

The procedure described above provides the weighted average and the covariance of the 
cross sections from different experiments contributing to a certain final state in a given 
range of energies. We now apply the trapezoidal rule. To perform the integration (5), we 
subdivide the integration range in fine energy steps and calculate for each of these steps 
the corresponding covariance (where additional correlations induced by the trapezoidal 
rule have to be taken into account). This procedure yields error envelopes between adja- 
cent measurements as depicted by the shaded bands in the corresponding figures. 

As a cross check, a different procedure of the evaluation of the integral has been ap- 
plied. For each final state, results of different experiments contributing to it in a given 
energy range are integrated separately using a rectangular method. After that a weighted 
average, based on the statistical and systematic errors combined in quadrature, is cal- 
culated. In some cases when correlations between systematic uncertainties of different 
experiments are known, they are taken into account after averaging the results with 
weights based on the statistical errors only. As mentioned above, if results of the inte- 
gration for different measurements are found to be inconsistent, the error is rescaled by a 



The difference between the results of the two described procedures is considered when 
estimating the systematic uncertainty on the numerical integration procedure. The sys- 
tematics also take into account variations of the energy interval where several data points 
are lumped into a single value, and the effect on the central value of the integral when 
including or not the correlations. The procedural systematics are added in quadrature to 
the experimental error on the integral. For instance, in the case of the 7r + 7r~ contribution, 
this procedural uncertainty amounts to 1.5 10~ 10 . 
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8 Specific Contributions 



In some energy regions where data information is scarce and reliable theoretical predictions 
are available, we use analytical contributions to extend the experimental integral. Also, 
the treatment of narrow resonances involves a specific procedure. 



8.1 The 7t + tt Threshold Region 

To overcome the lack of precise data at threshold energies and to benefit from the ana- 
lyticity property of the pion form factor, a third order expansion in s is used. The pion 
form factor F° is connected with the ir + n~ cross section via the expression 

KT = -^3 ^ • (28) 

The expansion for small s reads 

F% = 1 + l -{r 2 )^ s + Cl s 2 + c 2 s 3 + 0(s 4 ) . (29) 

Exploiting precise results from space-like data [128], the pion charge radius- squared is 
constrained to (r 2 ) n = (0.439 ± 0.008) fm 2 and the two parameters Ci j2 are fitted to the 
data in the range [2m n , 0.6 GeV]. In the case of r data, isospin corrections are taken into 
account as discussed before. 

The results of the fits are given in Table 3 and shown in Fig. 11. Good agreement 
is observed in the low energy region where the expansion should be reliable. Since the 
fits incorporate unquestionable constraints from first principles, we have chosen to use 
this parameterization for evaluating the integrals in the range up to 0.5 GeV. Systematic 
uncertainties due to the fitting procedure (fit boundaries, whether or not the coefficient 
c 2 is fixed) are small, albeit taken into account. 



8.2 Integration over the u and (j) Resonances 

In the regions around the u and resonances we have assumed in the preceding works that 
the cross section of the 7r + 7r~7r° production on the one hand, and the 7r + 7r~7r°, K + K~ as 
well as KgKl production on the other hand is saturated by the corresponding resonance 
production. In a data driven approach it is however more careful to directly integrate 
the measurement points without introducing prior assumptions on the underlying process 
dynamics [129]. Possible non-resonant contributions and interference effects are thus 
accounted for. 

Notwithstanding, a straightforward trapezoidal integration buries the danger of a bias: 
with insufficient scan density, the linear interpolation of the measurements leads to a 
significant overestimation of the integral when dealing with strongly concave functions 
such as the tails of Breit-Wigner resonance curves. This effect is particularly visible in 
the right hand plot of Fig. 12, showing the <f> resonance: the cross sections are measured 
by SND [46] (sum of the final states K + K~, KgK^ and -k + -k~it and corrected for missing 
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Data 


( ( \ o in r>i orrf" 


.Til 1 colli t 


Correlation matrix 


a had,LO ( 1Q -10) 














[2m v ± - 0.5 GeV] 




(r 2 ) 


CO 439 ± 008) fm 2 

\ KJ • -1- V t_y — 1 — KJ t KJ KJ KJ I 1111 


1 


* 








Cl 


(6.8 ± 1.9) GeV" 4 


-0.15 


1 




58.0 ± 1.7 ±l.l rad 




c 2 


(-0.7 ±6.8) GeV" 6 


0.09 


0.97 


1 








(0.439 ± 0.008) fm 2 


1 


* 






T 


Cl 


(3.3 ± 1.7) GeV" 4 


-0.15 


1 


* 


56.0 ± 1.6 ±0.3 SU (2) 




C2 


(13.2 ±5.7) GeV- 6 


0.09 


0.99 


1 




Table 3: 


Fit results 


of the low energy expansion (29) to 




and 


r data, the latter cor- 



rected for SU(2) breaking. The right column quotes the contributions to a^ ad ' LO ; integrated 
from threshold to 0.5 GeV. The errors are dominated by experiment, but take into account 
systematic uncertainties from the fitting procedure ( mainly the variation of the upper en- 
ergy cut yielding, e.g., an uncertainty of about 0.49 10~ 10 for r data). The systematics in 
a had,LO j rom radiative corrections (e + e~) and isospin breaking (r) (cf. Sections 4, 5.1) 
are quoted apart. 

modes, i.e., rescaled by (0.984 ±0.009)~ 1 [105]) and CMD-2 [48] {K Q S K° L only, rescaled by 
(0.337 ±0.005)~ 1 [105]). Shown in addition are the error band of the trapezoidal rule and 
the solution of a phenomenological fit of a BW resonance plus two Gaussians (only one 
Gaussian is necessary for the u, see left hand plot in Fig. 12) to account for contributions 
other than the single resonance. Both fits result in satisfactory x 2 values. Since we are 
only interested in the integral and do not want to extract dynamical parameters like 
phases or branching fractions, it is not necessary to parametrize the exact structure of 
the physical processes. We have accounted for the systematics due to the arbitrariness 
in the choice of the parametrization by varying the functions and parameters used. The 
resulting effects are numerically small compared to the experimental errors (see Table 4). 
It is clear from Fig. 12 that the fit function passes below the trapezoidal bands in the 
concave tails of both the and the u>. 

Table 4 gives the contributions to a^ ad > LO from the different energy domains covered by 
the experiments for both the uo and the 0. Since the experiments quote the cross section 
results without correcting for leptonic and hadronic vacuum polarization in the photon 
propagator (cf. the discussion in Section 4), we perform the correction here. Note that 
the data shown in Fig. 12 have been corrected for vacuum polarization. A small FSR 
correction (cf. Section 4) is applied to the results given in Table 4. The correction of 
hadronic vacuum polarization being iterative and thus only approximative, we assign half 
of the total vacuum polarization correction as generous systematic errors (cf. Section 4). 
In spite of that, the evaluation of a^ ad ' LO is dominated by the experimental uncertainties. 
Since the trapezoidal rule is biased, we choose the results based on the BW fits for the 
final evaluation of a^ ad ' L °. 
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s (GeV 2 ) 



Figure 11: Fit of the pion form factor from 4m 2 to 0.35 GeV 2 using a third-order Taylor 
expansion with the constraints at s = and the measured pion r.m.s. charge radius from 
space-like data [128]. The result of the fit is integrated only up to 0.25 GeV 2 . 

8.3 Narrow cc and bb Resonances 

The contributions from the narrow J/ip resonances are computed using a relativistic 
Breit-Wigner parametrization for their line shape. The physical values for the resonance 
parameters and their errors are taken from the latest compilation in Ref. [105]. Vac- 
uum polarization effects are already included in the quoted leptonic widths. The total 
parametrization errors are then calculated by Gaussian error propagation. This inte- 
gration procedure is not followed for the if)(3S) state which is already included in the R 
measurements, and for the T resonances which are represented in an average sense (global 
quark-hadron duality) by the bb QCD contribution, discussed next. 

8.4 QCD Prediction at High Energy 

Since the emphasis in this paper is on a complete and critical evaluation of spectral func- 
tions from low-energy data, we have adopted the conservative choice of using the QCD 
prediction only above an energy of 5 GeV. The details of the calculation can be found in 
our earlier publications [7, 11] and in the references therein. Only a very brief summary 
shall be given here. 

The perturbative QCD prediction uses a next-to-next-to- leading order 0{a\) expan- 
sion of the Adler D-function [130], with second-order quark mass corrections included [131]. 
R{s) is obtained by evaluating numerically a contour integral in the complex s plane. 
Nonperturbative effects are considered through the Operator Product Expansion, giv- 
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Figure 12: Cross sections of the uj (left) and <fi (right) resonances. The dots with error 
bars depict the measurements, the shaded band is the result of the trapezoidal rule within 
(correlated) errors and the function shows the phenomenological fit of a BW resonance plus 
one (uo) or two (cf)) Gaussians to account for other than the single resonance contributions. 
The bias of the trapezoidal rule when applied to a strongly concave (or convex) distribution 
is particularly visible in the tails of the (ft resonance when comparing to the BW fit. It 
leads to an overestimation of the integral. 



ing power corrections controlled by gluon and quark condensates. The value a s (M|) = 
0.1193±0.0026, used for the evaluation of the perturbative part, is taken as the average of 
the results from the analyses of r decays [6] and of the Z width in the global electroweak 
fit [132]. The two determinations have comparable uncertainties (mostly theoretical for 
the r and experimental for the Z) and agree well with each other. We conservatively take 
as final uncertainty the value quoted in either analysis. As for the other contributions, 
uncertainties are taken to be equal to half of the quark mass corrections and to the full 
nonperturbative contributions. 

A test of the QCD prediction can be performed in the energy range between 1.8 and 
3.7 GeV. The contribution to a^ ad > LO in this region is computed to be (33.87 ±0.46) 10~ 10 
using QCD, to be compared with the result, (34.9 ± 1.8) 10~ 10 from the data. The two 
values agree within the 5% accuracy of the measurements. 

In Ref. [11] the evaluation of a^ ad ' LO was shown to be improved by applying QCD 
sum rules. We do not consider this possibility in the present analysis for the following 
two reasons. First, it is clear that the main problem at energies below 2 GeV is now the 
inconsistency between the e + e~ and r input data, and this must be resolved with priority. 
Second, the improvement provided by the use of QCD sum rules results from a balance 
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_had,LO 

u n 
BW Fit 


(10 10 ) 
Trapez. 


a ( a had,LO) ( 1Q 

Exp. Fit 


10) 

BR 


Energy range 
(GeV) 


Type/Ref 




34.42 
2.51 


35.45 


0.63 
0.06 


0.37 
0.30 


0.27 
0.02 


0.760184 - 0.810 
0.300 - 0.760184 


CMD-2 [47] 
BWfit 


36.94 


37.96 


0.84 tot 


± 0.73 V p ± 0.30 FSR 


0.300 - 0.810 


Sum 


4> 


33.42 
32.84 


34.89 
34.28 


1.72 
0.72 


0.37 
0.39 


0.30 
0.49 


1.01017 - 1.03948 
1.01017 - 1.03948 


SND [46] 
CMD-2 [48] 


32.93 


34.37 




0.91tot 




1.01017 - 1.03948 


Average 


0.77 


1.10 


0.02 
0.06 


0.07 
0.01 


0.01 
0.01 


1 - 1.01017 
1.03948 - 1.055 


BW fit 
SND [46] 


34.80 


36.24 


0.92 tot 


±0.63 V p±0.14 FS R 


1 - 1.055 


Sum 



Table 4: Contributions to a^ ad > LO from the narrow resonances cj(782) (upper table) and 
0(1020) (lower table). Given are the results for the BW fit (first column) and the trape- 
zoidal rule (second column). The next three columns quote the experimental errors, the 
fit parameterization systematics and the uncertainty introduced by the correction for the 
missing decay modes of the resonances. The energy interval of the integration and the 
integration type (data or analytical function) are given in the last two columns. System- 
atic errors from the same sources, but for different energy regions are added linearly in 
the sum. All other errors are added in quadrature, the total errors being labelled 'tot 7 . 
Additional systematics are due to the vacuum polarization (VP) correction, taken to be 
half of the full correction, and to final state radiation (FSR) where the full correction is 
accounted as uncertain. 

between the experimental accuracy of the data and the theoretical uncertainties. The 
present precision of both e + e~ and r data, should they agree, is such that the gain would 
be smaller than before. This state of affairs will be reconsidered when the problems with 
the input data are sorted out. 

9 Results 

9.1 Lowest Order Hadronic Contributions 

Before adding up all the contributions to a^ ad ' LO , we shall summarize the procedure. On 
the one hand, the e + e~-based evaluation is done in three pieces: the sum of exclusive 
channels below 2 GeV, the R measurements in the 2-5 GeV range and the QCD predic- 
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tion for R above. Major contributions stem from the 2n (73%) and the two 4n (4.5%) 
channels. On the other hand, in the r-based evaluation, the latter three contributions are 
taken from r data up to 1.6 GeV and complemented by e + e~ data above, because the r 
spectral functions run out of precision near the kinematic limit of the r mass. Thus, for 
nearly 77% of a^ ad > LO (also contributing 80% of the total error-squared), two independent 
evaluations (e + e~ and r) are produced, the remainder being computed from e + e~ data 
and QCD alone. 

Fig. 13 gives a panoramic view of the e + e~ data in the relevant energy range. The 
shaded band below 2 GeV represents the sum of the exclusive channels considered in the 
analysis. It turns out to be smaller than our previous estimate [4], essentially because 
more complete data sets are used and new information on the dynamics could be incor- 
porated in the isospin constraints for the missing channels. It should be pointed out that 
the exclusive sum could lead to an underestimation of R, as some unmeasured higher 
multiplicity hadronic channels could start to play a role in the 2 GeV region. Neverthe- 
less, good agreement is observed at 2 GeV with the first inclusive data point from BES, 
thus indicating that the missing component is likely to be small. The QCD prediction is 
indicated by the cross-hatched band. It is used in this analysis only for energies above 
5 GeV. Note that the QCD band is plotted taking into account the thresholds for open 
flavour B states, in order to facilitate the comparison with the data in the continuum. 
However, for the evaluation of the integral, the bb threshold is taken at twice the pole 
mass of the b quark, so that the contribution includes the narrow T resonances, according 
to global quark-hadron duality. 

The contributions from the different processes in their indicated energy ranges are 
listed in Table 5. Wherever relevant, the two e + e~- and r-based evaluations are given. 
The discrepancies discussed above are now expressed directly in terms of a^ ad > LO giving 
smaller estimates for e + e~ data by (—21.2 ± 6.4 cxp ± 2.4 rad ± 2.6su(2) (±7.3 to tai)) 10~ 10 for 
the 2-7T channel and (—3.1 ± 2.6 cxp ± 0.3 rad ± 1.0 S u(2) (±2.9 to tai)) 10~ 10 for the sum of the 
4tt channels. The total discrepancy (-24.3 ± 6.9 cxp ± 2.7 rad ± 2.8 SU (2) (±7.9 to tai)) 10~ 10 
amounts to 3.1 standard deviations and precludes from performing a straightforward com- 
bination of the two evaluations. 



9.2 Results for a 



(i 



The results for the lowest order hadronic contribution are 



a 



had,LO 
had,LO 



(684.7 ± 6.0 cxp ± 3.6 rad ) 10 



-10 



a~~ = (709.0 ±5.1 cxp ±1.2 rad ±2.8 SU (2)) 10 



[e + e —based] 
10 [r-based] 



(30) 
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Figure 13: Compilation of the data contributing to a^ ad ' L °. Shown is the total hadronic 
over muonic cross section ratio R. The shaded band below 2 GeV represents the sum of the 
exclusive channels considered in this analysis, with the exception of the contributions from 
the narrow resonances which are given as dashed lines. All data points shown correspond 
to inclusive measurements. The cross-hatched band gives the prediction from (essentially) 
perturbative QCD, which is found to be in good agreement with the measurements in the 
continuum above 2 GeV. In this figure the bb threshold is indicated at the onset of BB 
states in order to facilitate the comparison with data in the continuum. In the actual 
calculation the threshold is taken at twice the pole mass of the b quark. 
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Modes 


Energy [GeV] 




a had,LO ( 10 "10) 
T (3) 


A(e+e" - r) 


LOW S exp. 7T + 7T~ 


[2m 7r ± — 


0.500] 


58.04 ± 1.70 ± 1.14 


56.03 ±1.61 ±0.28 


+2.0 ±2.6 


7T + 7r~ 


[0.500 - 


1.800] 


440.81 ± 4.65 ± 1.54 464.03 ± 3.19 ± 2.34 


-23.2 ±6.3 


7r°7, ?77 


[0.500 - 


1.800] 


0.93 ±0.15 ±0.01 




- 




[0.300 - 


0.810] 


36.94 ±0.84 ±0.80 




- 


7T + 7r~7r [below </>] 


[0.810 - 


1.000] 


4.20 ± 0.40 ± 0.05 




- 




[1.000 - 


1.055] 


34.80 ± 0.92 ± 0.64 




- 


7r + 7r 7T° [above </>] 


[1.055 - 


1.800] 


2.45 ± 0.26 ± 0.03 




- 


7T+7T-27r 


[1.020 - 


1.800] 


16.73 ±1.32 ±0.20 


21.44 ± 1.33 ±0.60 


-4.7 ± 1.8 


27T+27T" 


[0.800 - 


1.800] 


13.95 ±0.90 ±0.23 


12.34 ± 0.96 ± 0.40 


+1.6 ±2.0 


27r+27r-7r° 


[1.019 - 


1.800] 


2.09 ± 0.43 ± 0.04 


_ 


- 


^+^-37T° ( 2 ) 


[1.019- 


1.800] 


1.29 ±0.22 ±0.02 


_ 


- 


37r + 37r _ 


[1.350 - 


1.800] 


0.10 ±0.10 ±0.00 


_ 


- 


27r + 27T^27T 


[1.350 - 


1.800] 


1.41 ±0.30 ±0.03 


_ 


- 


tt+tt^tt ( 2 ) 


[1.350 - 


1.800] 


0.06 ± 0.06 ± 0.00 


_ 


- 


T/(— > 7T + 7T~7, 27)7T + 7T~ 


[1.075 - 


1.800] 


0.54 ±0.07 ±0.01 


_ 


- 


7r°7)^° 


[0.975 - 


1.800] 


0.63 ±0.10 ±0.01 


_ 


- 


7r°7)(vr^) 


[1.340 - 


1.800] 


0.08 ± 0.01 ± 0.00 


_ 


- 




[1.055 - 


1.800] 


4.63 ± 0.40 ± 0.06 


- 


- 




[1.097- 


1.800] 


0.94 ±0.10 ±0.01 




- 


tfOtfi^T (2) 


[1.340 - 


1.800] 


1.84 ±0.24 ±0.02 


- 


- 


ifi?7T° (2) 


H A An 
[1.44U — 


i cnnl 
l.oUUJ 


0.60 ±0.20 ±0.01 




- 


KKtttt ( 2 ) 


[1.441 - 


1.800] 


2.22 ±1.02 ±0.03 




- 


R = J2 excl. modes 


[1.800 - 


2.000] 


8.20 ±0.66 ±0.10 




- 


it [Data] 


[2.000 - 


3.700] 


26.70 ± 1.70 ±0.00 






J/Y> 


[3.088 - 


3.106] 


5.94 ±0.35 ±0.03 




- 


^(25) 


[3.658 - 


3.714] 


1.50 ±0.14 ±0.00 






i? [Data] 


[3.700 - 


5.000] 


7.22 ±0.28 ±0.00 






Rudsc [QCD] 


[5.000 - 


9.300] 


6.87 ±0.10 ±0.00 






Rudscb [QCD] 


[9.300 - 


12.00] 


1.21 ±0.05 ±0.00 






Rudscbt [QCD] 


[12.0 - 


- oo] 


1.80 ±0.01 ±0.00 






J] (e + e _ — > hadrons) 


[2m 7r ± 


— oo] 


684.7 ± 6.0 exp 
± 3.6 rad 


709.0±5.1 cxp 

± l-2 ra d ± 2.8 SU ( 2 ) 


-24.3 ± 7.9tot 



1 Not including u and <fi resonances (see text). 

2 Using isospin relations (see text). 

3 e + e~ data are used above 1.6 GeV (see text). 



Table 5: Summary of the ajj ' contributions from e + e~ annihilation and r decays. The 
uncertainties on the vacuum polarization and FSR corrections are given as second errors 
in the individual e + e~ contributions, while those from isospin breaking are similarly given 
for the r contributions. These 'theoretical' uncertainties are correlated among all channels, 
except in the case of isospin breaking which shows little correlation between the 2n and 
47r channels. The errors given for the sums in the last line are from the experiment, the 
missing radiative corrections in e + e~ and, in addition for r , SU(2) breaking. 
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Figure 14: Comparison of the results (31) with the BNL measurement [16]. Also shown 
are our previous estimates [3, 11] obtained before the CMD-2 data were available, and the 
recent evaluation of Hagiwara et al. [129]. 



Adding the QED, higher-order hadronic, light-by-light scattering and weak contributions 
as given in Section 2, the results for obtained 



7 SM 

V 



11 659 169.3 ± 7.0 had ± 3.5 LBL ± 0.4 qed+ ew) 10~ 10 
11 659 193.6 ± 5.9 had ± 3.5 LBL ± 0.4 QED+EW ) KT 10 



e e 



-based] 
[r— based] . 



These values can be compared to the present experimental average given in Eq. (2). 
Adding experimental and theoretical errors in quadrature, the differences between mea- 
sured and computed values are found to be: 



exp 



„exp n SM 



(33.7 ± 11.2) 1(T 10 
(9.4 ± 10.5) 10" 10 



e 1 e 



-based] 
[t— based] , 



(31) 



corresponding to 3.0 and 0.9 standard deviations, respectively. A graphical comparison 
of the results (31) with the experimental value is given in Fig. 14. Also shown are our 
previous estimates [3, 11] obtained before the CMD-2 and the new r data were available 
(see discussion below), and the recent evaluation of Hagiwara et al. [129]. 
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10 Discussion 



10.1 The Problem of the 2tt Contribution 

The significant discrepancy between the e + e~ and r evaluations of a^ ad ' LO is a matter 
of concern. In this section we comment on the relevant aspects of the problem. Since 
our earlier work [4, 7, 11] was based on a combined analysis of e + e~ and r data, we feel 
important to summarize the main changes (all expressed in 10~ 10 units) in the dominant 
2n contribution where the r contribution makes its impact: 

• the new CMD-2 data [22] produce a downward shift of the e + e~ evaluation by 1.9 
(well within errors from previous experiments), while the final error is reduced from 
±12.5 to ±5.1 with an additional ±2.4 from missing radiative corrections, 

• the new ALEPH data [24] increases the r evaluation by 3.5, which is within the 
previous experimental uncertainty of ±7.2 (in our previous analyses, we did not 
quote the results of a r-based analysis alone, but only those from the combined 
spectral functions), 

• including the CLEO data in the r evaluation improves the precision, but further 
raises the central value by 4.0, 

• although including the OPAL data has little effect on the overall precision, it also 
increases the result by 1.9, 

• the new complete isospin symmetry-breaking correction, including the re-evaluation 
of the Sew factor, increases the r evaluation by 0.2 with respect to the previous 
one [4, 7, 11], which is well within the quoted error of ±2.5. 

The previous (unpublished) difference between the e + e~- and r-based evaluations of the 
27r contributions, A (a^ ad ' LO ) = — 10.8±12.5exp,ee±7.2 cxPiT ±2.5su(2), was consistent 

with zero, allowing the two spectral functions to be combined into an improved common 
estimate. 

In spite of the fact that every change was within its previously estimated errors, the 
two results are not consistent anymore so that one must address the question of the 
possible origin of the problem. 

In principle, the observed discrepancy for the 2ir contribution, (— 21.2±7.3), or (— 4.2± 
1.4)% when expressed with respect to e + e~, could be caused by any (or the combination 
of several) of the following three effects which we examine in turn: 

• The normalization of e + e~ data 

Here, as below, 'normalization' does not necessarily mean an overall factor, but refers 
to the absolute scale of the 'bare' cross section at each energy point. There is no 
cross check of this at the precision of the new CMD-2 analysis. The only test we can 
provide is to compute the e + e~ integral using the experiments separately. Because of 
the limited energy range where the major experiments overlap, we choose to perform 
the integration in the range of y/s from 610.5 to 820 MeV. The corresponding 
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contributions are: 313.5 ± 3.1 for CMD-2, 321.8 ± 13.9 for OLYA, 320.8 ± 12.6 
for CMD, and 323.9 ± 2.1 for the isospin-corrected r data. No errors on radiative 
corrections and isospin breaking are included in the above results. 

• The normalization of r data 

The situation is quite similar, as the evaluation is dominated by the ALEPH data. 
It is also possible to compare the results provided by each experiment separately, 
with the spectral functions normalized to the respective hadronic branching ratios. 
Leaving aside the region below 500 MeV where a fit combining analyticity con- 
straints is used, the contributions are: 460.1 ±4.4 for ALEPH, 464.7 ±9.3 for CLEO 
and 464.2 ±8.1 for OPAL, where the common error on isospin breaking has been 
left out. The three values are consistent with each other and even the less precise 
values are not in good agreement with the e + e~ estimate in this range, 440.8 ± 4.7, 
not including the error on missing radiative corrections. This is in line with the con- 
clusion drawn from the comparison of branching ratios presented in Fig. 10. The 
larger values obtained with the CLEO and OPAL spectral functions are related to 
their relatively higher level below the p resonance, as can be observed in Fig. 3. 

At the level of the r — > z/ r 7r _ 7r branching ratio, which controls the normalization 
of the 7r~7r° spectral function, stringent tests can be applied to the ALEPH results. 
We stress the fact that the branching fractions are obtained by a global procedure 
where all r decay final states are considered, down to branching ratios of a few 10 -4 , 
from a very clean initial sample [24]. The most critical part in the analysis is the 
separation of channels with different ir° multiplicities. The it~it final state could 
be spoiled from the adjacent channels n~ and 7r~27r° by inadequate understanding 
of the 7 identification and the 7r° reconstruction. The observed branching ratios 
for these two modes are in agreement with expectations, based for the first one 
only on the assumption of universality of the /i — r couplings in the weak charged 
current (which is tested at the 3 10~ 3 level using the r electronic branching ratio 
and the lifetime), and for the second one on the isospin relation with the 2-k~it + 
branching ratio: B n - = (-0.08 ± 0.11 cxp ± 0.04 th )%, and B n27T o - B^ = 
(±0.06 ± 0.17 cxp ± 0.07 t h)%- These two tests provide confidence that the precise 
determination of the branching ratio for r — > u t tv~tv is on solid ground, as the 
observed discrepancy would require a shift of 1.1% on this quantity. 

Apart from an overall normalization effect, differences could originate from the shape 
of the measured spectral functions. If all three spectral functions are normalized 
to the world average branching ratio (our final procedure), then the results for the 
contribution above 0.5 GeV become: 459.9 ±3.6 for ALEPH, 465.4 ±5.1 for CLEO 
and 464.5 ±5.1 for OPAL, with a common error of ±2.4 from the 7T7r and leptonic 
branching ratios and the uncertainty on isospin breaking left out. Again the results 
are consistent and their respective experimental errors give a better feeling of the 
relative impact of the measurements. 

• The isospin-breaking correction applied to r data 

The basic components entering SU(2) breaking have been identified. The weak 
points before were the poor knowledge of the long-distance radiative corrections and 
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the quantitative effect of loops. Both points have been addressed by the analysis of 
Ref. [30] showing that the effects are small and covered by the errors previously ap- 
plied. The overall effect of the isospin-breaking corrections (including FSR) applied 
to the 2ir r data, expressed in relative terms, is (—1.8 ±0.5)%. Its largest contribu- 
tion (—2.3%) stems from the uncontroversial short-distance electroweak correction. 
Additional contributions must be identified to bridge the observed difference. 

One could question the validity of the chiral model used. The authors of Ref. [30] 
argue that the corrections are insensitive to the details of their model and essentially 
depend only on the shape of the pion form factor. As the latter is known from 
experiment to adequate accuracy, there seems to be little space for improvement. 

Thus we are unable at this point to identify the source of the discrepancy. More ex- 
perimental and theoretical work is needed. On the experimental side, additional data is 
available from CMD-2, but not yet published. As an alternative, a promising approach 
using e + e~ annihilation events with initial state radiation (ISR), as proposed in Ref. [133], 
allows a single experiment to cover simultaneously a broad energy range. Two experimen- 
tal programs are underway at Frascati with KLOE [134] and at SLAC with BABAR [135]. 
The expected statistics are abundant, but it will be a challenge to reduce the systematic 
uncertainty at the level necessary to probe the CMD-2 results. However, the experimental 
technique being so different, it will be in any case valuable to compare the results with 
the present ones. As for r's, the attention is now focused on the forthcoming results from 
the B factories. Again, the quality of the analysis will be determined by the capability to 
control systematics rather than the already sufficient statistical accuracy. On the theory 
side, the computation of more precise and more complete radiative corrections both for 
e + e~ cross sections and r decays should be actively pursued. 

10.2 Other Points of Discussion 

Other points are worth to be discussed: the 4n spectral functions, the ui and <f> resonances 
and the sum of exclusive channels from 1.6 to 2 GeV. 

As already pointed out in Sections 6.1 and 6.2, the quality of both e + e~ and r data in 
the 2n + 2'K~ and 7r + 7r~27r° final states is not as good as for the 7r + 7r~ channel. Agreement 
is observed in the former channel at the 10% level, while in the latter a large discrepancy 
is found (see the values in Table 5) which to this level cannot be attributed to isospin 
breaking. Since significant differences are found within the e + e~ data sets, we feel that it 
is a priority to clarify the experimental situation in this sector. The ISR program being 
conducted with the BABAR experiment should be able to shed some light upon this 
problem [135]. 

Compared to previous estimates, the uj and <p resonance contribution is now directly 
evaluated with the measured cross section, rather than integrating a Breit-Wigner function 
computed with averaged parameters. The uo value is basically unchanged, while a large 
downward shift of 4.3 has been found for the <fi contribution. The origin of this change 
lies in the fact that the recent measurement of the <p lineshape yields a total width which 
is significantly smaller: decreased by Qa in the last two years [105]! 
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Revisiting the situation of the exclusive channels in the 1.6-2 GeV range has led to 
significant changes, the origin of which are twofold: (i) more information (obtained in the 
study of r decays, see Section 3.1) on the decay dynamics in the 6n channel could be used 
to bound the 7r + 7r~47r contribution, and (n) some data with poor quality were discarded, 
resulting in smaller contributions in the 37r + 3"7r~ and 27r + 27r~27r° channels. As a result, 
and unlike the conclusion reached in Refs. [136, 129], we find the sum of the exclusive 
processes to be in reasonable agreement with the inclusive measurements of R in this 
range [86] . At any rate it is clear that better data should be taken in this energy region. 
The BABAR ISR physics program should be able to make an important contribution here 
as well. 

Due to the last two points which are only relevant to e + e~ data — the contribution 
from the <fi resonance and the multi-pion channels — our new evaluation comes out to be 
significantly smaller than before. 

10.3 Comparison to Other Evaluations of a^ ad ' L0 

Here we restrict our discussion to recent evaluations which have been published since 1998, 
i.e. Refs. [137, 138]. Previous estimates were considered in our earlier publication [11]. 
A common feature of Refs. [137, 138] is that they use both e + e~ and r data for the 2n 
contribution. However, their analyses are based on the preliminary CMD-2 data [43] which 
are not corrected for vacuum polarization and FSR. Because of this, they fail to notice 
the discrepancy between e + e~ and r data. In addition, no mention of isospin symmetry 
breaking, and how to correct for it, is made in Ref. [137], shedding some doubts about 
the validity of the combination of the two data sets. The relatively high values obtained 
in these analyses compared to the present one are due in part to these problems. 

The recent analysis of Hagiwara et al. [129] does include the final CMD-2 data. Our 
e + e^-based result agrees with their evaluation using inclusive hadron production for en- 
ergies above 1.6 GeV. However, as pointed out before, our re-evaluated sum of exclusive 
channels in this range is consistent within errors with the inclusive rate. 

10.4 Consequences for a(M^) 

In spite of the fact that the present analysis was focused on the theoretical prediction 
for the muon magnetic anomaly, it is possible to draw some conclusions relevant to the 
evaluation of the hadronic vacuum polarization correction to the fine structure constant 
at Mf . The problem found in the 2ir spectral function is less important for Aa(Mf) 
with respect to the total uncertainty, because the integral involved gives less weight to 
the low-energy region. The difference between the evaluations using the 2tt, 4tc and 27r27r° 
spectral functions from e + e~ and r data are found to be: 

AcCi(^I) " A< d (Mf ) = (-2.79 ± 0.43 cc ± 0.26 rad ± 0.55 r ± 0.30 su(2) ) 1(T 4 . 

S v ' 

0.80 

While this low-energy contribution shows a 3.5 standard deviation discrepancy (when 
adding the different errors in quadrature), it also exceeds the total uncertainty of 1.6 1CT 4 
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on Aa(Mf) which was quoted in Ref. [11]. It is worth pointing out that such a shift 
produces a noticeable effect for the determination of the Higgs boson mass M H in the 
global electroweak fit [139]. With the present input for the electroweak observables [132] 
from LEP, SLC and FNAL yielding central values for Mh around 100 GeV, going from the 
e + e~ to the r-based evaluation induces a decrease of Mh by 16 GeV using all observables 
and by 20 GeV when only the most sensitive observable, (sin 2 #w)eff, is used. 

11 Conclusions 

A new analysis of the lowest-order hadronic vacuum polarization contribution to the 
muon anomalous magnetic moment has been performed. It is based on the most recent 
high-precision experimental data from e + e~ annihilation and r decays in the 7T7T channel. 
Special attention was given to the problem of isospin symmetry breaking and the cor- 
responding corrections to be applied to r data. A new theoretical analysis of radiative 
corrections in r decays was used and found to be in agreement with previous estimates. 
A complete re-evaluation of the contributions of e + e _ annihilation cross sections in the 
energy range up to 2 GeV has been performed. Incorporating the recently corrected contri- 
bution from light-by- light scattering diagrams, the full prediction for the muon magnetic 
anomaly is obtained. 

The main results of our analysis are the following: 

• the new evaluation based solely on e + e~ data is significantly lower than previous es- 
timates and is in conflict with the experimental determination of by 3.0 standard 
deviations. 

• the new precise evaluations of the dominant nn contributions from e + e~ annihilation 
and isospin-breaking corrected r decays are not anymore in agreement with each 
other. A discussion has been presented for possible sources of the discrepancy which 
could not be resolved. This situation is a matter of great concern, as the r-based 
prediction of is in better agreement with the experimental value, from which it 
deviates by non-significant 0.9 standard deviations. 

More experimental and theoretical work is needed to lift the present uncertainty on 
whether or not new physics has been uncovered with the muon magnetic moment. 
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